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ABSTRACT 

Context. At the Very Large Telescope Interferometer, the purpose of the fringe-tracker FINITO is to stabilize the optical path differences 
between the beams, allowing longer integration times on the scientific instruments AMBER and MIDI. 

Aims. Our goal is to demonstrate the potential of FINITO for providing //-band interferometric visibilities, simultaneously and in addition to 
its normal fringe-tracking role. 

Methods. We use data obtained during the commissioning of the Reflective Memory Network Recorder at the Paranal observatory. This device 
has permitted the first recording of all relevant real-time data needed for a proper data-reduction. 

Results. We show that post-processing the FINITO data allows valuable scientific visibilities to be measured. Over the several hours of our 
engineering experiment, the intrinsic transfer function is stable at the level of ±2%. Such stability would lead to robust measurements of 
science stars even without the observation of a calibration star within a short period of time. We briefly discuss the current limitations and the 
potential improvements. 

Key words. Techniques: high angular resolution - Techniques: interferometric - Instrumentation: high angular resolution - Instrumentation: 
interferometers 



1. Introduction 

The purpose of the FINITO facility at the Very Large Telescope 
Interferometer (VLTI, see Scholler et al. 2006) is to compen- 
sate for the random optical delay (OPD) introduced by atmo- 
spheric turbulence. If not corrected, such perturbation would 
make the interferometric fringes jitter on the detector, pre- 
venting the use of long integration time on scientific instru- 
ments. FINITO was offered to the astronomical commu nity in 



Septe mber 2007, and since then several AMBER (Petr ov et al 



2007h programs have already benefite d from the enhanced data 
quahty produced by fringe-tracking ( Le Bouquin et al. I l2008bt 
Wittkowski etai] |2008). 

Nevertheless, FINITO can potentially tackle another im- 
portant issue at VLTI: measuring robust, absolutely cal- 
ibrated, interferometric visibilities in the near-IR. Indeed, 
AMBER has been designed to be a highly 



rewarding 



spectro-interferometer (Petrov et al. 2007; Weigelt et al. 2007; 
[Pomici a no de Souza et a l. 2007; Meilland et al...2007.) . but this 
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has required technical trade-offs. It typically reads out more 
slowly than an atmospheric coherence time, and so is critically 
sensitive to fringe stabilization, and hence to seeing and me- 
chanical vibrations. This issue is exacerbated by VLTI service- 
mode only providing a few calibration stars per night, which 
is insufficient for properly sampling the time variations of the 
instrumental response (called transfer function). AH together, 
this mak es the absolute cahb ration of AMBER fringe contrast 
difficult dMillour et alj|2007h . 

We have recently carried out the necessary developments 
to make FINITO able to provide visibility measurements. No 
hardware modifications were needed as FINITO already inter- 
ferometrically combined the light of 2 or 3 telescopes to stabi- 
lize the OPD. The main issue was the high frame rate of the sys- 
tem running generally at more than IkHz. This was overcome 
by the commissioning of the so-called VLTI Reflective Memory 
Network Recorder (RMNrec), which listens to the communica- 
tions between several systems of VLTI and is able to store the 
real-time data into proper FITS files. 

In this paper, we investigate the possibility of post- 
processing these real-time data to measure interferometric vis- 
ibilities and to properly calibrate them. In Sect. |2] we summa- 
rize the FINITO and RMNrec instrumental setup, the test ob- 
servations made, and the simple data reduction used for this 
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Fig. 1. Sketch of the Finito optical bench showing the fibered 
Fringe Modulator Units (FMU) and the two co-axial beams 
combiners. Black solid lines are for the mirrors, while dashed 
lines are for semi-transparent devices, {^nj 
interferometric outputs and {P\,P(),P2] are the photometric out- 
puts. The real-time flux at those seven outputs are recorded in 
real time in a camera not displayed here. 



demonstration. In Sect. [3] as a classical test, we compare the 
measured angular diameter of some well-known stars with pre- 
viously published values. We discuss the results, the accuracy 
obtained, and the perspectives for reaching higher precision in 
the future. The paper ends with brief conclusions and perspec- 
tives. 

2. Observations and data reduction 

2.1. RMNRec and FINITO instrumental setup 

For our study, we use the commissioning data from the 
RMNrec machine at Paranal. By looking at the Reflective 
Memory Network at high frequency, this computer constructs 
FITS files with the real-time data of several VLTI guiding and 
tracking facilities: error vector of the centroids as measured on 
the IRIS guiding camera; real-time fringes, phase and flux as 
measured on FlNlTO; and fringe-tracking loop status (search, 
idle, lock) and OPD offsets as handled by the OPD-Controller 
software (OPDC). Our long-term goal is to benefit for all in- 
formation in the data processing of VUTl instruments. For the 
purpose of this paper, we only focus on the FINITO and OPDC 
data because they are of scientific interest in themselves, since 
they contain the interferometric visibility of the observed tar- 
get. 

A de tailed des c riptio n of the FINITO facility can be 



found in iGai et al.l (120041) . Recent updates can be found in 



Le Bouquin et al.l ( 2008 aV Brieflv. FINITO interferometric ally 
combines two pairs of telescopes from a telescope triangle in 
the //-band (see Fig. [U. The optical path length of two of the 
three beams are modulated internally using a piezo-driven ex- 
pansion ring changing the length of the fibers, called a Fringe 
Modulator Unit (FMU). We generally use a triangular modula- 



Target 


HD 


H 


Spectral 


Published in H 


Name 


Number 


(mag) 


Type 


(mas) 


a Hya 


81797 


-0.99 


K3II 


9.44 ± ? 


N Vel 


82668 


-0.35 


K5III 


6.89 ± 0.08 


5 Vir 


112300 


-1.0 


M3III 


9.91 ±0.1 


BK Vir 


108849 


-0.37 


M7III 


10.73 ± 0.23 


a Hya 


73840 


2.0 


K3III 


2.62 ±0.1 


^ Vir 


110014 


2.17 


K2III 


1.95 ± 0.02 


i Hya 


83618 


0.95 


K2.5III 


3.39 ± 0.05 


eCrt 


99167 


1.19 


K5III 


3.45 ± 0.04 



Note: Spectral types and H-band magnitudes have been extracted 
from the Simbad a nd 2MASS database, and diameters from 
CHARM2 catalog jRichichi et alj2005l) . 



tion of 5 //-band fringes in amplitude peak-to-peak and a sam- 
pling of 4 readouts per fringes. Each modulation ramp is here- 
after called a scan. Photometric signals are extracted by Glan- 
Taylor polarizers, while the other polarization component is fed 
to a three-way, pairwise, co-axial beam combiner. The central 
beam (labelled 0) is split into equal intensity fractions and then 
combined with the modulated beams (respectively 1 and 2). 
This produces two pairs of complementary interferometric out- 
puts (tt phase difference), for the combination 0-1 and 0-2, re- 
spectively. The 3 photometric signals (P\,Pq,P2) and the 4 in- 
terferometric signals (^fo, ^10,^20^20) routed to a Rockwell 
PICNIC camera, which is read at a frequency between 500Hz 
and 2kHz. These output signals are processed in real time to 
measure the fringe phases, used by the OPDC to reject fringe 
motion for frequencies below 20Hz. 

2.2. Observations 

Commissioning RMNrec data were taken during the nights of 
2008-01-28, 2008-01-29, and 2008-01-30 with the interfero- 
metric baseline AO-DO alone. FINITO was therefore tracking 
on a single channel (beams 0-2, known to provide better per- 
formances). Sequences of 3 files of 30s were recorded at each 
telescope pointing. Over the 3 nights, this leads to a total of 1 .1 
hours of fringe integration over approximately 5 hours of ob- 
servation (about 45 pointing). Neither the AMBER nor MIDI 
instruments were used simultaneously, so the only interfero- 
metric data available are from FINITO. 

The seeing ranged between 0.75" and 2.2" during the obser- 
vations. The fringe-tracking loop was always able to lock the 
fringes with average performances: locking time greater than 
80% and residual OPD generally less than 0.3 A^m RMS. The 
fringes were never lost for more than a second. The main dif- 
ference between good and bad files are the number of fringe 
jumps from which the tracking loop had to recover 

Concerning the sources, well known small and large stars 
were observed sequentially to test the feasibility of astrophys- 
ical data processing by measuring stellar diameters. In the fol- 
lowing we respectively call the small and large stars calibrators 
and targets. Their stellar parameters are summarized in Table[Tl 
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Fig. 2. Example of a FINITO real-time fringe signal C2o(f) (upper panel) on a short sequence of 0.2s, with its associated FMU 
modulation (red curve in lower panel). The FMU modulation was 5 fringes in amplitude, with 4 readouts per fringes, and 0.5ms 
per readout, leading to 10ms per scan. The scan numbering is given at the bottom of each plot. Gray fringes in the lower panel are 
the theoretical //-band fringe-packet of scans 5 and 17, assuming a infinite linear modulation ramp. The part of the packet that 
is effectively modulated is over-plotted in black. This pattern is replicated for each scan (mirrored in time when the FMU slope 
is negative). After the fringe jump (f - 22.68s), the fringe-packet appears shifted with respect to the middle of the modulation 
window. This explains the apparent oscillation in the envelope of signal C2o(f)- 



2.3. Data reduction 

Step 1 : The FINITO raw data are processed in real time in 
order to compute the correction send to the Delay Lines. In this 
process, the two raw interferometric outputs are corrected and 
calibrated from the photometric fluctuations by the following 
classical formula (here for beams 2-0): 

/i> - K2^2()+ Pi - ^0^20+^0 



r+ - 



( 



'20 



2 y/K2^2Q+ kq^2Q+ Pi Po 
- Ki^nr Pi - kq^iq- Po 



20 



(2) 



^20- ^0^20- Pi Pa 

The /c-matrix is the internal transmission ratio between the pho- 
tometric and the interferometric outputs. It is calibrated each 
time we change star by measuring the flux at the outputs while 
observing one beam at a time. 

The two cleaned signals are then subtracted to enhance the 
signal-to-noise ratio (SNR): 

C20+ - C20- 



C20 



(3) 



By doing so, all the correlated noises between the two outputs 
should disappear, while the fringes, perfectly opposed in phase, 
are amplified. The resulting signal C2o(f) is processed in real 
time to compute the phase residual and the group-delay, which 
is used to feed the fringe-tracking loop. It is accessible on the 
RMN and can be stored into FITS files by the RMNrec ma- 
chine. An example of fringe sequence is displayed in Fig. |2] 
For illustration purpose, we chose a sequence with an obvious 
fringe jump (at t = 22.68s). In the last 7 scans, the fringe-packet 
was not centered inside the modulation window. Those fringe 
jumps are caused by temporary flux dropouts and can occur at 
unpredictable times. Between a fringe jump and the time the 
OPDC compensates for it, the average fringe contrast is obvi- 
ously reduced. All these scans should be efficiently detected 
and discarded before computing the average visibility. 



Step 2: We used two simple criteria for the scan selection: first 
the fringe-tracking loop should be closed during the complete 
scan, and secondly the fringe-packet maximum should be less 
than one fringe from the scan center. Additionally, the scans 
were weighted according to their average SNR. To do this, 
we configured the RMNrec to store all required OPDC signals 
(loop status, fringe-packet position, SNR). 

We then computed the power spectral density (PSD) of 
each of the accepted scans, averaged all of them, and integrated 
the power around the expected fringe frequency: 



^700Hz 

I^L. = < \Ca(t)\^ > df - bias 

J300HZ 



(4) 



where <> denotes the average over the accepted scans. This 
power is a quadratic estimation of the fringe visibility, assum- 
ing the bias is properly removed. In our data set, we estimated 
it by interpolating the values at the edges of the fringe peak; see 
Fig. [3]for typical PSDs with their bias removed. We decided to 
integrate the fringe power in the frequency range 300-700Hz 
to take into account the small amount of energy spread due to 
the shape of the fringe-packet in the direct space (broad-band 
fringes). We checked that integrating over a wider range 200- 
800Hz does not significantly change the final calibrated visibil- 
ities. 

All reduced visibility points firaw are displayed in Fig. H] 
On the largest stars, the super-synthesis eff'ect is clearly visi- 
ble, making the baseline lengths (and therefore the visibilities) 
slowly change over time. However, these raw visibilities can- 
not be used directly for estimating the stellar diameters. They 
have to be calibrated from the interferometric response of the 
whole instrumental chain, also called the transfer function. 
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Fig. 3. Average power spectral density of typical FINITO se- 
quences of 30s, after scan selection and weighting (see Step 3). 
The square visibility is obtained by integrating the fringe power 
around 500Hz (filled region), after removing the bias (dashed 
red curve). 
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Fig. 4. Raw visibiUties for all observations, along with the date 
(month-day) and the detector integration time (italic) informa- 
tion. See Fig.|6]for the color code of each source. 



Step 3: We computed the instantaneous transfer function yuo 
by dividing the observed visibilities measured on each calibra- 
tors with those expected according to the diameters of Table [T] 



(5) 



In this well-known formula, B is the baseline length of each vis- 
ibility point fi, A the eff'ective wavelength, and the Airy func- 
tion (Fourier transform of a uniform disk). In the following, we 
assume an effective wavelength of 1.625yL(m, as measured by 
Bauvir et al. (private communication). Results are plotted in 
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Fig. 5. Top: Transfer function computed for each calibrator ob- 
servation, along with the date (month-day) and the detector in- 
tegration time (italic) information. Middle and bottom: same 
transfer function points but plotted versus the FINITO resid- 
ual phase RMS (accepted scans only) and versus the Paranal 
DIMM seeing. See Fig.|6]for the color code of each source. 



Fig. |5] Given the overall stability of the transfer function, we 
decided to average all points of a given setup, defined as the 
same night and the same detector integration time (DIT). We 
used this average transfer function to calibrate all observed vis- 
ibiUties of the corresponding setup: 



(6) 



<;"0 > 



The resulting calibrated visibilities are plotted against the base- 
line length in Fig. |6] To extract uniform disk diameters for 
each target, we fit the full //(B) data set with the expression 
yi(0B//i). Measured uniform disk diameters are summarized 
in Table |2] as well as comparison of previously published val- 
ues. 
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Fig. 6. Calibrated visibilities of all observations plotted ver- 
sus the baseline length. The solid lines are the literature di- 
ameters used to compute the transfer function (for the calibra- 
tors), or the best-fitted diameters of the data set (for the targets). 
The shaded regions show the diameter uncertainties. Symbols 
are for the different setups, defined as same night and same 
FINITO integration time (DIT). 

Table 2. Published versus measured diameters. 



Target 


Published 


Method 


Lbd 


Measured 


Nb. of 


Name 


(mas) 


(') 


(jim) 


in H (mas) 


pointing 


a Hya 


9.09 ± 0.09 


LBI 


0.8 


8.95 ±0.15 


13 


N Vel 


6.89 ± 0.08 


IND 


1.6 


6.6 ± ? 


1 


5 Vir 


9.91 ±0.10 


LBI 


0.8 


10.1 ±0.1 


4 


BKVir 


10.73 ±0.23 


LBI 


2.2 


10.2 ±0.1 


5 



Indirect (IND) or long baseline interferometry (LBI) estimation. 



affecting all scans in the same ways, or bias in the averaging 
process. We see at least four potential origins: 



1. In long baseline interferometry from the ground, the 
transfer function depends on the atmospheric turbulence 
strength. To explain the observed variations, we should in- 
voke fast changes in the atmospheric turbulence at the level 
of one minute or less. But we do not see any correlations 
when plotting the transfer function estimations against the 
seeing or the fringe-jitter (see Fig.|5]), the latter being a good 
proxy for the coherence time. Moreover, we expect FINITO 
to be robust against turbulence fluctuations because it com- 
bines the use of single-mode fibers, high sampling rate, and 
real-time piston correction. 

2. The /f-matrix and its associated eiTors are not stored in the 
FITS files produced by RMNrec, and are not propagated 
through Eqs.[T]to|6] Because this /c-matrix is only recom- 
puted at each pointing, any errors on it would lead to sys- 
tematic offsets between sequence of 3 consecutive files on 
a given star. If present, such behaviors are not dominant in 
Fig.|5] demonstrating that errors on the /c-matrix alone can- 
not explain the observed discrepancies. 

3. As explained in Step 2 of Sect. 12.31 we discarded the scans 
where the central fringe was not centered in the modulated 
window. The efficiency of this selection process is corre- 
lated to the accuracy on the fringe-packet center, which 
may depend on the observing conditions. 

=A-. As seen in the PSD plotted in Fig. [3] the fringe 
power around 500Hz is probably contaminated by a low- 

_ frequency signal. Its exact spectral shape is unknown and 
looks unstable, making it difficult to fit correctly. In the 
data reduction, we used a linear interpolation between ar- 
bitrary frequencies to estimate and remove it (see Step 2 of 
Sect. I2.3l l. This may have introduced an additional error in 
the integrated fringe power. 



3. Results and discussion 

All four measured diameters are within 5% of the pre- 
viously published diameters. This was the goal of the 
study: the good agreement proves the capability of the 
FINITO/OPDC/RMNrec facility to provide meaningful inter- 
ferometric visibility measurements. Interestingly, the transfer 
function is stable over several hours. It is one of the most im- 
portant parameters in obtaining robust visibility measurements. 

3.1. Visibility accuracy 

Total uncertainty on our visibility points is difficult to estimate. 
We used the bootstrapping method to compute the contribu- 
tion of the statistical dispersion inside each file. This technique 
has the important advantage of not assuming a Gaussian dis- 
tribution of the measurements. Resulting statistical errors are 
very small, on the order of ±0.1%, twhich is inside the sym- 
bol size of Fig. |5] This does not explain the observed disper- 
sion between consecutive data points. The total uncertainty is 
therefore dominated by systematic effects, that is to say, errors 



3.2. Perspectives toward liiglier accuracy 

At the time of our observations, RMNrec was only able to 
store the signal C2o{t), not the FINITO raw signal themselves. 
Therefore, as explained in Sect. 12.31 the results presented in 
this paper are based on data that were previously processed in 
real time by Eqs.[T]and[3] Data reduction would highly bene- 
fit from a logging of the raw FINITO data. For instance, pho- 
tometric calibration during post-processing can make use of 
Wiener Filtering, which is not possible in real time since this 
filter is not causal. Additionally, logging the /c-matrix and asso- 
ciated errors is mandatory to check its stability and its effective 
contribution into the final error 

FINITO is a fringe tracker, meaning that fringes are locked 
in phase. Using coherent average over several scans (instead of 
power-average as in this paper) can reduce the impact of the 
bias on the final result. On the other hand, the transfer func- 
tion of coherent-average depends on the fringe-tracking loop 
performances. All together, this approach is worth looking at 
since it may be a useful trade-off for fainter stars, where bias 
subtraction becomes more complicated. We plan to do so as 
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soon as the logging of the OPDC/FINITO data becomes a stan- 
dard procedure. 

3.3. Perspectives toward fainter targets 

More generally, we should discuss the expected performances 
on fainter targets. Our main concern is related to Step 2 of 
the proposed data-reduction (discarding the fringe packets that 
were not properly centered inside the modulated window). 
Indeed, the efficiency of such scan selection, while good on 
high-SNR fringes, most probably degrades rapidly on noisy 
data. The simplest and most robust solution is to modulate 
an OPD-range significantly wider than the frin ge-packet size. 



as with the instrument VINCI for instance dKervella et al 



2004 |Le Bouqui n et al.1 120041) . However, the FINITO fringe- 
tracking capability would be lost because the phase measure- 
ments become extremely noisy while outside the fringe-packet. 
Performing packet-tracking would still be possible, but this 
prohibits the use of long integration on the scientific camera 
(required for medium and high spectral resolution modes of 
AMBER). Because of its potential scientific interest on faint 
targets, we are looking forward to implementing long scans and 
group-tracking as an additional mode of FINITO. 

4. Conclusions 

We have demonstrated the possibility of post-processing the 
real-time fringe-tracking data of the Very Large Telescope 
Interferometer for astronomical purposes. We proposed a sim- 
ple strategy for reducing the data, based on scan selection and 
power integration. We were able to recover the //-band diame- 
ters of four already-known giant stars. The good general agree- 
ment proved the capability of the FINITO/OPDC/RMNrec fa- 
cility to provide meaningful interferometric visibility measure- 
ments. 

Interestingly, the transfer function level is dominated by 
the integration time. It is stable over several hours, with pos- 
sible significant changes from one day to the next. We have not 
been able to find any obvious correlation with the atmospheric 
seeing and fringe-tracking performances, even if we agree that 
our dataset does not allow us to draw definite conclusions. 
Post-processing the FINITO/OPDC data therefore provides a 
rough, but robust, estimation of the //-band absolute visibility 
even without the observation of a calibration star close in time. 
This complements the scientific observations very well, gener- 
ally done in another spectral range (near-IR for AMBER, and 
A^-band for MIDI) and with higher spectral resolution (up to 
10000 for AMBER), but with only poor absolute calibration 
of the continuum visibility. We emphasize that this additional, 
crucial measure is obtained simultaneously and without any ad- 
ditional overhead time. 

This demonstration is also promising in the context of pre- 
cision interferometry at the VLTI. We were able to reach accu- 
racy of +2% on the raw visibilities with a simple data reduction 
method. We showed that the overall accuracy is most proba- 
bly limited by systematic eiTors. Statistical errors only about 
+0. 1 % for 30s integration on H = 2™ unresolved targets with 
the 1.8m Auxiliary Telescopes. 



Based on the success of this demonstration, we now plan 
to routinely log and archive the VLTI fringe-tracking data. In 
this context, we emphasize the need for a dedicated public data 
reduction pipeline. 
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